Uncertainties of clock and shift operators for an electron in one-dimensional nonuniform lattice systems
Gong Long-Yan1, 2, 3, †, Ding You-Gen2, Deng Yong-Qiang2
Department of Applied Physics, Nanjing University of Posts and Telecommunications, Nanjing 210003, China
Institute of Signal Processing and Transmission, Nanjing University of Posts and Telecommunications, Nanjing 210003, China
National Laboratory of Solid State Microstructures, Nanjing University, Nanjing 210093, China

 

† Corresponding author. E-mail: lygong@njupt.edu.cn

Project supported by the National Natural Science Foundation of China (Grant Nos. 61475075 and 61170321).

Abstract

The clock operator U and shift operator V are higher-dimensional Pauli operators. Just recently, tighter uncertainty relations with respect to U and V were derived, and we apply them to study the electron localization properties in several typical one-dimensional nonuniform lattice systems. We find that uncertainties ΔU2 are less than, equal to, and greater than uncertainties ΔV2 for extended, critical, and localized states, respectively. The lower bound LB of the uncertainty relation is relatively large for extended states and small for localized states. Therefore, in combination with traditional quantities, for instance inverse participation ratio, these quantities can be as novel indexes to reflect Anderson localization.

1. Introduction

The Heisenberg uncertainty principle is one of the most important ideals in quantum mechanics.[1] It is often expressed by no-go statements.[2] In a general expression, two noncommuting observables do not simultaneously have well-defined values. This property has no classical analogue, so it is a basic difference between classical and quantum mechanics. As is well known, the principle plays a fundamental role in quantum theory, such as nonlocality,[3] quantum measurements,[4] and quantum phase transitions.[5] Furthermore, it has important applications in current quantum technology, for instance, entanglement detection,[6] quantum cryptography,[7] and quantum metrology.[8]

In fact, the principle was first proposed by Heisenberg,[1] and originally formulated precisely by Kennard[9] as ΔxΔp/2, where the uncertainties are quantified by standard deviations of the position and momentum. Subsequently, Robertson[10] generalized it to arbitrary pairs of noncommuting observables A and B, which is given by ΔAΔB ≥ 〈Ψ|[A,B]|Ψ〉/2 in state |Ψ〉. Just recently, several uncertainty relations with respect to the clock operator U and shift operator V were derived.[11] The expressions of U and V are given by[1114] where [Z] is the integer part of Z. Here, and are complete orthonormal bases in d-dimensional Hilbert space, respectively. The two bases are related by the discrete Fourier transform (DFT) Based on the DFT, the operators U and V can be rewritten as and . In fact, the expressions of U and V in Eqs. (B2) and (B4) of Ref. [11] are originally from Eqs. (5) and (7) of Ref. [14]. They should be the same. However, they are not completely the same. We have checked them and found that there exist typo errors in Eqs. (B2) and (B4) of Ref. [11]. Physically, the operators U and V are the non-Hermitian generalization of the traditional 2-dimensional Pauli operators to higher-dimensional ones.[12,13] The definitions of clock operator U and shift operator V are independent of the Hamiltonian, so it can be applied to general quantum systems.

The uncertainties of U and V are respectively given by[11,14] where |Ψ〉 is the wave function of a quantum state and 0 ≤ ΔU2 ≤ 1 (0 ≤ ΔV2 ≤ 1). One of the tighter uncertainty relations that they obey is[11] where Δ(3) = 〈Ψ|ΨU〉〈ΨU|ΨV〉〈ΨV|Ψ〉, |ΨU〉 = U|Ψ〉, |ΨV〉 = V|Ψ〉 and Φ = arg Δ(3). We define the right side in the inequation is the lower bound

On the other hand, Massar and Spindel pointed out that the uncertainties ΔU2 and ΔV2 can measure the localization properties of quantum states.[14] However, they do not give the corresponding explicit relation. Besides, Bagchi and Pati pointed out that LB depends on quantum states.[11] However, they did not study the relation between the values of LB and state localization properties. In this paper, we fill the gaps just mentioned. In the Anderson localization phenomenon,[1520] a state |Ψ〉 may be localized, critical or delocalized, which depend on models and model parameters. There are many typical and extensively studied one-dimensional models, such as the slowly varying potential model,[2123] the deterministic nearest-neighbor tight-binding models with exact mobility edges,[24] and the random-dimer potential model.[2527] It is interesting how these quantities, including ΔU2 and ΔV2 and LB, are related to Anderson localization.

The rest of the paper is organized as follows. In Section 2, we introduce the tight-binding Hamiltonian for a single electron in one-dimensional (1D) lattice systems. In Section 3, we analyze the behaviors of ΔU2 and ΔV2, and LB in three typical models. Finally, we give a summary of our main results.

2. The tight-binding Hamiltonian for a single electron in one-dimensional lattice systems

A single electron in 1D lattice systems is considered. The corresponding tight-binding Hamiltonian is described by where εj is the on-site potential on the jth site, t is a nearest-neighbor hopping integral, and d is the lattice size. Further, , where is the creation operator of jth site and |0〉 is the vacuum. Models are specified by the choice rules of {εj}. The system can also be described by the tight-binding equation The general wave function, i.e., eigenstate |Ψβ〉 with eigenenergy Eβ for Hamiltonian in Eq. (7) can be written as the superposition where is the amplitude of the βth wave function at jth site.

3. Numerical results

In numerical calculations, for different models, we directly diagonalize the tight-binding equation (8) with the periodic boundary condition at finite system sizes and obtain all eigenenergies Eβ and corresponding eigenstates |Ψβ〉. From Eqs. (3), (4), and (6), we obtain ΔU2 and ΔV2, as well as LB. For simplicity and convenience, based on Eqs. (3) and (4), we define a reduced relative uncertainty

3.1. Slowly varying potential model

For the slowly varying potential (SVP) model,[2123] the on-site potential in Eq. (7) can be written as Here, λ is the potential strength and α is a real number; φ is the original phase, which is redundant in localization. It is found that for 0 < υ < 1, there are two mobility edges (MEs) at Ec = ± (2.0 − λ)t provided that λ < 2.0t. Extended states are in the middle of the band (|Eβ| < |Ec|). Localized states are at the band edges (|Eβ| > |Ec|); states are more localized as eigenenergies Eβ are nearer band edges. At MEs, states are critical.

We take λ = 1.0t as an example. At the same time, πα = 0.2 and υ = 0.7 are set, which are the same as that in Refs. [21]–[23]. In this case, two MEs are Ec = ±t. We denote and , respectively. States are extended for |Eβ| < t and localized for |Eβ| > t. Figure 1(a) shows that at extended states, uncertainties ΔU2s are near ones. At localized states, ΔU2s decrease as states vary from MEs to band edges; at the same eigenenergies Eβ, the larger lattice sizes, the smaller ΔU2s are. There are sharp transitions in ΔU2s at MEs, where states are critical. Thus ΔU2s can well characterize the localization properties of the three different kinds of states in position space. At the same time, figure 1(b) shows that on the whole, uncertainties ΔV2s are decreasing as states vary from band center to band edges. There is a deep pit in ΔV2s when eigenenergies Eβ are near each ME. The values of ΔV2 are almost independent of the lattice sizes Ns. Fortunately, figure 1(c) shows that the reduced relative uncertainty Ur > 0 for extended states and Ur < 0 for localized ones. For larger lattice size N, Ur almost equals to zero for critical states (seen in inset). It means that uncertainties ΔU2s are greater than, equal to, and less than uncertainties ΔV2s for extended, critical, and localized states, respectively. Therefore, different kinds of states can be distinguished from each other with the comparisons of ΔU2 and ΔV2.

Fig. 1. (color online) For the slowly varying potential model, (a) ΔU2, (b) ΔV2, (c) Ur, (d) LB, and (e) IPR as functions of eigenenergies Eβ, respectively. Here, λ = 1.0t, πα = 0.2, and υ = 0.7; the vertical dashed lines are for the functions and , respectively; the horizontal dashed line in panel (c) is for Ur = 0; lattice sizes d = 1000, 2000, and 10000, respectively. Inset in panel (c): partial enlarger for Eβ near .

Figure 1(d) shows that for extended states, LBs are near ones. For localized states, on the whole, the more localized the states are, the smaller the LBs are. There are rapid transitions in LBs as states move from the extended into the localized regime. A deep pit is in LBs when eigenenergies Eβ are near each ME, so LBs depend on states, and different kinds of states can also be distinguished from each other based on LB.

Furthermore, compared to our results intuitively, we also calculate the inverse participation ratio , which is a traditional and often used quantity to characterize localization properties of states. Roughly, it is proportional to the localization length of states.[28] Figure 1(e) shows that for extended states, IPRs are relatively large. For localized states, IPRs are relatively small, and at a fixed Eβ, the larger the d is, the smaller IPR is. On the whole, all our studied quantities are in good correspondence with IPR.

3.2. Deterministic nearest-neighbor tight-binding models with exact mobility edges

One family of the models are specified by on-site potentials[24] The on-site potential is a smooth function of in the open interval (−1,1). Each value of corresponds to a different tight-binding model. is the original phase, which is redundant in localization. For a quasiperiodic modulation, b is set to be irrational. There exist self-duality points under a generalized duality transformation.[24] Based on it, the analytical critical condition about MEs when is where sgn(x) is the sign function.

In Fig. 2, we take the parameters of λ = −0.9t and as examples. From Eq. (15), an ME Ec = t/3. In this case, all states are localized when eigenenergies Eβ > Ec, and almost all states, except some states at subband edges, are extended when Eβ < Ec. In calculations, the irrational number b in Eq. (14) is approximated by the ratio of successive Fibonacci numbers: Fj = Fj − 1 + Fj − 2, with F0 = F1 = 1. In this way, choosing b = Fj − 1/Fj and lattice size d = Fj, we can obtain the periodic approximant for the quasiperiodic potential.

Fig. 2. (color online) For the deterministic nearest neighbor tight binding models with exact mobility edges, (a) ΔU2, (b) ΔV2, (c) Ur, (d) LB, and (e) IPR as functions of eigenenergies Eβ, respectively. Here, λ = −0.9t and ; the vertical dashed line is for the function Eβ = Ec = t/3; the horizontal dashed line in panel (c) is for Ur = 0; lattice sizes d = 987, 2584, and 10946, respectively.

Similar to that for the SVP model, figure 2 shows that for extended states, uncertainties ΔU2s are near ones, the reduced relative uncertainty Ur > 0, and LBs are near ones. For localized states, ΔU2s are near zeros, Ur < 0, LBs are relatively small. There are good correspondences for them with IPRs. So all these quantities can characterize the localization properties of different kinds of states in the model.

3.3. Random-dimer potential model

In the random-dimer potential (RDP) model,[2527] the on-site potential in Eq. (7) is given by where Va and Vb are dimer energies. The potential is random but with short correlations, which is different from potentials in the above two models where potentials are determined (without disorder). It was found[25] that delocalized (extended) states exist when |VaVb| < 2.0t and there are only localized states when |VaVb| > 2.0t. At finite lattice sizes d, the number of delocalized states is roughly proportional to ; for infinite lattice sizes, there are only states with eigenenergy Eβ = Va(Vb) are delocalized.[26]

In Fig. 3, we take Va = −Vb = −0.75t as examples. At the same time, we choose q = 0.5, which corresponds to the most random situation. In this case, states with eigenenergies Eβ near Va(Vb) are delocalized. The results of Fig. 3 are obtained for average of 500, 250, and 50 samples for lattice sizes d = 1000, 2000, and 10000, respectively. The averages with more samples give the same results.

Fig. 3. (color online) For the random-dimer potential model, (a) ΔU2, (b) ΔV2, (c) Ur, (d) LB, and (e) IPR as functions of eigenenergies Eβ, respectively. Here, Va = −Vb = −0.75t; the vertical dashed lines are for the functions Eβ = Va and Eβ = Vb, respectively; the horizontal dashed line in panel (c) is for Ur = 0; lattice sizes d = 1000, 2000, and 10000, respectively.

Similar to that for the SVP model and the deterministic nearest-neighbor tight-binding models with exact mobility edges, figure 3 shows that for states with eigenenergies Eβ very near Va(Vb) (delocalized states), ΔU2s are near ones, Ur > 0, LBs are also near ones. For other states (localized ones), at the same Eβ, ΔU2s decease with lattice sizes Ns; Ur < 0 and LBs are relatively small. There are good correspondences between these quantities and IPRs. Thus, they can reflect the localization properties of different kinds of states in the model.

4. Discussion and conclusions

It is known that the operators U and V are the non-Hermitian generalization of the Pauli operators from two dimensions to higher dimensions. The simplest higher-dimensional quantum states may be states for an electron in lattice systems. In the three studied models, we found that ΔU2 can well characterize state localization properties. The reduced relative uncertainty Ur can also do so. At the same time, Resta proposed and used a U-related quantity to study many-body localization.[29] Here, we found that U itself can well reflect the localization properties of single electron lattice systems. Besides, Bagchi and Pati pointed out that LB [Eq. (6)] depends on states.[11] Here, we give further results that LBs are relatively large for extended states and relatively small for localized states. Therefore, all these studies present more pictures than existing works.

For the three studied models, we found that uncertainties ΔU2 > ΔV2, ΔU2 = ΔV2, and ΔU2 < Δ V2 for extended, critical, and localized states, respectively. A similar result has been found by contrast of Shannon information entropies of states for an electron in position and momentum space,[30] or by contrast of the distinguishability and visibility for an electron in lattice systems.[31] Thus, these quantities can be used as novel indexes to characterize Anderson localization.

Reference
[1] Heisenberg W 1927 Z. Phys. 43 172
[2] Busch P Heinonen T Lahti P 2007 Phys. Rep. 452 155
[3] Oppenheim J Wehner S 2010 Science 330 1072
[4] Deutsch D 1983 Phys. Rev. Lett. 50 631
[5] Sachdev S 1999 Quantum Phase Transitions Cambridge Cambridge University Press
[6] Horodecki R Horodecki P Horodecki M Horodecki K 2009 Rev. Mod. Phys. 81 865
[7] Scarani V Bechmann-Pasquinucci H Cerf N J Dušek M Lütkenhaus N Peev M 2009 Rev. Mod. Phys. 81 1301
[8] Giovannetti V Lloyd S Maccone L 2011 Nat. Photon. 5 222
[9] Kennard E H 1927 Z. Phys. 44 326
[10] Robertson H P 1929 Phys. Rev. 34 163
[11] Bagchi S Pati A K 2016 Phys. Rev. 94 042104
[12] Schwinger J 1960 Proc. Natl. Acad. Sci. USA 46 570
[13] Sylvester J J 1909 The Collected Mathematics Papers of James Joseph Sylvester Cambridge Cambridge University Press III
[14] Massar S Spindel P 2008 Phys. Rev. Lett. 100 190401
[15] Anderson P W 1958 Phys. Rev. 109 1492
[16] Abrahams E 2010 50 Years of Anderson Localization Singapore World Scientific
[17] Zhu H J Xiong S J 2010 Chin. Phys. 19 037107
[18] Gong L Y Tong P Q 2008 Chin. Phys. 17 0674
[19] Liao J Shi G Liu N Li Y Q 2016 Chin. Phys. 25 117201
[20] Chou C Ho C 2014 Chin. Phys. 23 110302
[21] Das Sarma S He S Xie X C 1988 Phys. Rev. Lett. 61 2144
[22] Griniasty M Fishman S 1988 Phys. Rev. Lett. 60 1334
[23] Thouless D J 1988 Phys. Rev. Lett. 61 2141
[24] Ganeshan S Pixley J H Das Sarma S 2015 Phys. Rev. Lett. 114 146601
[25] Dunlap D H Wu H Phillips P W 1990 Phys. Rev. Lett. 65 88
[26] Gangopadhyay S Sen A K 1992 J. Phys.: Condens. Matter 4 9939
[27] Datta P K Giri D Kundu K 1993 Phys. Rev. 47 10727
[28] Kramer B Machinnon A 1993 Rep. Prog. Phys. 56 1469
[29] Resta R 1998 Phys. Rev. Lett. 80 1800
[30] Gong L Wei L Zhao S Cheng W 2012 Phys. Rev. 86 061122
[31] Gong L Xue B Li W Cheng W Zhao S 2016 Phys. Rev. 94 032125